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摘 要 : 针对 现 有 医学 图 像 加 密 算 法 在 加 密 效率 和 安全 性 上 的 不 足 ， 提 出 一 种 基于 2D sine Logistic 混沌 映射 的 医学 
图 像 小 波 域 加 密 算 法 。 算 法 首先 利用 整数 小 波 变换 将 医学 图 像 从 空域 转换 为 频 域 ， 充 分 打破 像素 间 的 相关 性 ; 其 次 ， 
利用 2D sine Logistic 混沌 映射 生成 混沌 序列 ,选取 三 级 小 波 分 级 的 低频 系数 LL3 进行 扩散 和 置 乱 加 密 ， 提 高 加 密 效 
和: 天 且 半 级 小 分 的 中 光 系 铸 本 2 和 LH2 过 入 认 加 密 ， 有 加密 因 人 让 在 的 明基 站 是; 后 
加 密 后 的 小 波 系数 进行 小 波 逆 变换 得 到 加 密 图 像 。 实 验 仿真 结果 表明 ， 算 法 具有 高 安全 性 和 加 密 效率 ， 与 现 有 空 
方法 相 比 ， 加 密 时 间 约 为 1/40; 与 现 有 频 域 方法 相 比 ， 在 保证 加 密 效 率 情 况 下 具有 更 好 的 加 密 图 像 隐蔽 性 
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Abstract: Aiming at the shortage of medical images encryption algorithm in efficiency and security, this paper proposed a 


novel medical image encryption algorithm in frequent domain based on 2D sine Logistic chaotic mapping. Firstly, the 
proposed algorithm used the integer wavelet transform of LeGall5/3 to transform medical image from frequent domain to 
spatial domain. In this way, the high correlation among pixels has disappeared. Then, the proposed algorithm utilized 2D 


和 sine logistic chaotic mapping to generate the chaotic Sequences. In order to improve the encryption efficiency, this algorithm 


chose the low frequency coefficients of three-level decomposed to encrypt with chaotic diffusion and scrambling 
mechanism. In the same time, this algorithm also selected the medium-high frequency coefficients of two-level decomposed 
to encrypt by chaotic diffusion mechanism for solving the problem of obvious outline in encrypted image. Finally, this 
algorithm transformed the encrypted wavelet coefficients by the inverse integer wavelet transform and obtained the 
encrypted image. Experimental simulation results demonstrate that the proposed algorithm is able to protect medical images 
with low time complexity and a high security level. Compared with the existing spatial method, it reduces the encryption 
time by about 0.025 times. In addition, the proposed algorithm has better image imperceptibility as well as an excellent 
encryption efficiency compared to the existing frequent method . 
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被 研究 者 们 证 明 具 有 高 的 安全 性 和 加 密 效率 广 9。 

前 ， 基 于 混沌 映射 的 医学 图 像 加 密 算法 可 以 分 为 两 大 
随 着 “互联 网 + 医疗 ”的 不 断 发 展 , 医学 图 像 的 在 线 传输 变 ”类 : 基于 空域 和 基于 频 域 的 方法 ， 前 者 利用 混沌 映射 生成 的 
得 越 来 越 普遍 。 数 字 医 学 图 像 在 临床 医疗 诊断 中 起 着 至 关 重 ”混沌 序列 对 医学 图 像 的 像素 进行 加 密 ， 典 型 的 方法 如 文献 
要 的 作用 ,但 由 于 其 含有 病人 的 重要 隐私 信息 ， 其 安全 问题 [7~9]。 文 献 四 设计 了 一 个 基于 混沌 映射 的 置 乱 与 替换 扩散 架 
也 受到 研究 者 们 的 广泛 关注 03。 在 众多 的 安全 保护 方法 中 ， ” 构 用 于 医学 图 像 加 密 ， 对 所 有 像素 进行 异 或 加 密 后 改变 其 位 
加 密 仍 然 是 对 数据 进行 主动 保护 的 有 效 方法 。 相 比 自然 图 像 ， ” 置 ， 算 法 具有 高 的 加 密 性 能 ， 但 由 于 医学 图 像 的 数据 量 大 ， 
医学 图 像 具 有 数据 量 大 、 典 型 的 分 区 域 特征 、 高 的 像素 相关 ”算法 的 加 密 效 率 不 高 。 为 了 解决 基于 空域 算法 的 加 密 效 率 问 
性 、 直 方 图 分 布 不 均匀 等 特点 ， 一 些 适 用 于 自然 图 像 的 加 窗 。 题 ，Moumen 等 人 局 提 出 选取 图 像 的 部 分 像素 进行 加 密 ， 利 
方法 并 不 适应 于 医学 图 像 ， 如 传统 的 ASE、3DES 等 算法 并 用 图 着 色 理 论 选 取 图 像 中 的 部 分 像素 进行 混沌 加 密 。Pareek 
不 能 满足 大 数据 量 的 医学 图 像 实时 加 密 。 混 沌 映射 具有 伪 随 ” 等 个 根 据 医学 图 像 的 分 区 域 特性 ， 提 出 对 医学 图 像 的 感 兴 
机 性 、 人 遍历 性 和 初 值 敏感 性 等 特点 ， 根 据 混沌 映射 产生 的 混 芒 区 域 (前 景区 域 》 进 行 混沌 加 密 。 但 上 述 方法 需要 特征 抽 
沌 序列 具有 安全 密 钥 的 良好 特性 ， 混 沌 密码 学 逐渐 成 为 密码 取 、 模 式 匹 配 和 其 他 先 验 知识 ， 降 低 了 算法 的 可 操作 性 ， 如 
学 一 个 新 的 研究 方向 ， 利 用 混沌 映射 进行 医学 图 像 加 密 已 经 。” 感 兴趣 区 域 提取 时 需要 对 区 域 进行 分 割 、 由 于 区 域 的 不 规则 
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特性 需要 对 区 域 进行 表示 L9H。 基于 频 域 的 加 密 方法 首先 将 医 。 较 差 ， 如 Logistic 混沌 系统 已 经 证 明 存 在 着 安全 缺陷 叫 。 高 
学 图 像 从 空域 转换 为 频 域 ， 然 后 对 选取 的 频 域 系数 进行 加 密 。 ” 维 混沌 系统 具有 更 为 复杂 的 系统 结构 和 混沌 性 能 ， 但 实现 费 
后 进行 逆 变 换 得 到 加 密 图 像 ， 典 型 的 方法 如 文献 [11~13]。 文 ”时 。Hua 等 人 [5 将 两 个 一 维 的 混沌 系统 Logistic 和 Sine 映 系 
献 [11] 提 出 了 一 种 基于 余弦 变换 和 混沌 映射 的 加 密 方 法 ， 结合 ， 提 出 了 2D sine Logistic 混沌 映射 ， 并 且 证 明了 其 在 安 
有 较 高 的 加 密 效率 ， 但 由 于 余弦 变换 涉及 到 浮 点 数 运算 ， 在 ”全 性 和 实现 效率 上 相 比 其 他 高 维 混沌 系统 具有 优势 。2D sine 
不 增加 专门 浮 点 运算 硬件 的 情况 下 保证 医学 图 像 的 加 解密 可 ”Logistic 混沌 映射 可 表示 为 

道 性 是 不 现实 的 。 为 了 解决 这 一 问题 ，Wu 等 人 02 和 深 洗 青 I @) 
等 人 031 均 提出 了 基于 整数 小 波 变换 和 超 混 沌 映射 的 加 密 方 yi = 0(sin(Axin) + PB)y,(1— y,) 

法 ， 整 数 小 波 变换 避免 了 系数 运算 中 的 浮 点 数 问题 ， 而 超 混 其 中 :xy es[0，we[0] ，Be[0,3] ， 当 ws[0.905.1 并 且 8 接 
沌 映射 提升 了 算法 的 安全 性 ， 算 法 仅 对 小 波 分 解 后 的 低频 系 ” 近 3 时 ，2D sine logistic 混沌 映射 具有 超 混沌 行为 。 

数 进行 混沌 加 密 ， 大 大 提高 了 加 密 效 率 ， 但 超 混沌 序列 的 生 | 2 “所 
成 在 时 间 效 率 上 是 瓶颈 。 由 于 频 域 系数 与 图 像 空域 像素 之 间 ey 


bz 


没有 明确 的 对 应 关系 ， 频 域 系数 的 改变 往往 会 带 来 大 量 空域 得 直方 向 中 蝇 关 了 第 
像素 值 的 变化 ， 选 择 频 域 系数 进行 图 像 加 密 具 有 较 好 的 安全 Lm | am 


特性 ， 在 加 密 效率 上 高 于 空域 算法 。 

综 上 所 述 ， 基 于 空域 的 加 密 算法 相 比 频 域 方法 ， 具 有 更 
高 的 加 密 性 能 ， 但 加 密 效 率 较 低 ;基于 频 域 的 方法 仅 对 部 分 ee 
系数 进行 加 密 ， 为 了 提高 算法 的 安全 性 ， 需 要 设计 更 加 复杂 


的 混沌 映射 ,但 超 混沌 系统 的 设计 对 算法 的 加 密 效率 是 挑战 ， 

并 且 通 过 对 文献 [13] 进 行 实验 分 析 发 现 ， 仅 对 低频 系数 进行 人 | 

加 密 ， 在 医学 图 像 存在 明显 边界 的 情况 下 加 密 图 像 中 轮廓 显 到 1 图 像 三 级 小 波 分 级 的 系数 子 带 及 分 解 实 例 

现 比 较 严 重 ， 算 法 的 安全 性 降低 。 针 对 上 述 问 题 ， 本 文 提出 Fig. 1 Coefficient subband of three-level wavelet transform and 
基于 2D sine Logistic 混沌 映射 的 医学 图 像 频 域 加 密 算 法 ， 算 decomposition example of medical image 

法 采用 的 2D sine Logistic 混沌 映射 已 经 被 证 明 具 有 超 混沌 系 1.3 问题 提出 

统 特性 ， 并 且 由 简单 的 低 维 混沌 系统 组 合 而 成 ， 具 有 更 加 简 文献 [15] 首 先 利 用 2D sine Logistic 混沌 映射 生成 混沌 序 
单 的 系统 结构 。 另 外 ， 对 小 波 分 解 后 的 中 高 频 子 带 进行 简单 ” 列 ， 然 后 利用 混沌 序列 设计 像素 位 置 置 乱 和 扩散 算法 ， 算 法 
的 扩散 加 密 ， 打 破 原 有 图 像 中 边缘 和 纹理 信息 ， 解 决 加 密 图 接 对 图 像 的 所 有 像素 进行 操作 。 表 1 给 出 了 文献 [15] 方 法 


y gy 的 结果 可 以 看 言 息 入 接近 理想 值 8， 说 明了 

口 与 问题 担 
1 ”相关 知识 与 问题 提出 法 良好 的 加 密 效 果 ， 但 随 着 图 像 尺寸 的 增 大 ， 加 密 所 需 的 
1.1 整数 小 波 变换 时 间 大 幅 增 加 。 虽 然 更 优 的 实验 环境 可 能 带 来 更 好 的 实验 结 
为 了 保证 图 像 加 解密 的 可 逆 性 ， 采 用 基于 提升 方案 的 整 ” 果 ， 但 对 于 大 数据 量 的 医学 图 像 来 说 ， 算 法 并 不 适用 。 一 方 


像 中 的 轮廓 问题 。 在 自然 图 像 Lena 不同 的 图 像 尺 寸 ) 中 的 测试 结果 ， 从 表 1 
中 
算 


数 小 波 变换 LeGall5/3 ， 该 变换 将 整数 映射 到 整数 ， 并 且 小 波 面 由 于 算法 在 空域 进行 ， 需 要 加 密 的 数据 量 大 ， 男 外 一 方面 
的 正 变换 和 逆 变 换 过 程 中 保证 了 信息 的 无 损 性 ， 是 JPEG200 于 设计 的 扩散 和 置 乱 算法 过 于 复杂 。 
压缩 标准 中 指定 的 可 逆 小 波 变换 0 。 设 一 个 一 维 的 信和 号 表 1 文献 [15] 方 法 在 Lena 中 的 测试 结果 
s=[s01,502,…,Sow] ， 其 中 50 为 整数 ,可 采用 式 (1) 的 小 波 正 变 换 Table ] Lena’s testing results in reference 15 
将 信号 进行 一 级 小 波 分 级 。 载体 图 像 加 密 图 像 信息 粹 加密 时 间 /s 解密 时 间 /s 
六 -| 2(502, + so222) +112| Lena(128*128) 7.9882 0.610 0.594 
re- ey +|1/ 4 4d +112| (1) Lena(256*256) 7.9972 8.14 8.09 
Lena(S512*512) 7.9993 40.044 41.538 
其 中 : 4 为 高 频 系 数 ，s5, 为 低频 系数 ， 对 于 一 级 小 波 分 级 得 图 2 给 出 了 载体 图 像 ct 采 用 文献 [13] 方 法 加 密 后 的 结果 。 
到 的 低 高 频 系 数 可 以 进一步 利用 式 (1) 进 行 更 高 级 的 分 解 ， 如 ”文献 [13] 方 法 仅 对 三 级 小 波 分 解 的 低频 系数 加 密 ， 从 图 1 中 
5 代表 三 级 小 波 分 级 的 低频 系数 。 对 于 二 维 的 图 像 和 矩阵, 可 ”可 以 看 出 ， 得 到 的 图 像 出 现 明 显 的 轮廓 。 事 实 上 ， 通 过 实验 
以 把 图 像 的 行 和 列 向 量 分 别 进 行 小 波 分 解 ， 对 于 分 解 后 的 系 ” 分析， 对 于 轮廓 对 比 不 分 明 的 医学 图 像 ， 文 献 [13] 具 有 较 好 
数 算 当 


E 阵 ， 可 采用 式 (2) 的 小 波 道 变换 得 到 原始 图 像 。 的 加 密 效果 ， 当 医学 图 像 中 前 景 和 背景 区 域 对 比 度 较 大 ， 呈 
人 现 黑 白 分 明 时 加 密 图 像 出 现 明显 轮廓 。 加 密 图 像 中 一 旦 存在 
a (2) ”原始 图 像 的 轮廓 ， 会 严重 影响 加 密 算法 的 安全 ， 攻 击 者 可 以 

通过 选择 明文 攻击 来 破译 加 密 算法 。 


故 力 : j 为 分 解 级 数 。 图 1 给 出 了 图 像 三 级 小 波 分 解 后 的 频 2 ”算法 描述 
域 系数 子 带 和 医学 图 像 分 解 实例 。 每 级 小 波 分 解 得 到 四 个 频 ee 


率 子 带 ,分 别 为 LL、HL、LH 和 HH， 其 中 LL 为 低频 子 带 ， ”2.1 算法 模型 

集中 了 原始 图 像 的 绝 大 部 分 能 量 ， 是 原 图 的 逼近 子 图 ，HL 本 文 算法 模型 如 图 3 所 示 。 算 法 首先 将 原始 医学 图 像 进 

和 LH 分 别 为 垂直 方向 和 水 平方 向 的 中 高 频 子 带 ， 代 表 原 图 行 三 级 小 波 正 变 换 ， 得 到 其 系数 矩阵 ， 其 中 选取 三 级 分 解 低 

的 纹理 和 边缘 信息 ， 为 图 像 轮廓 信息 的 主要 集中 地 。 频 系 数 LL3 和 二 级 分 解 中 高 频 系 数 LH 和 HL 作为 待 加 密 系 

1.2 2D sine Logistic 混沌 映射 数 , 其 他 系数 保持 不 变 ; 然后 在 给 定 混 沌 初 值 和 参数 情况 下 ， 
维 的 混沌 系统 具有 更 为 简单 的 系统 结构 ， 但 安全 性 能 “利用 2D sine logistic 混沌 映射 生成 混沌 序列 1 和 2; 接着 将 
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系数 LL 利用 序列 1 进行 扩散 加 密 ， 并 将 加 密 后 的 结果 利用 
序列 2 进行 置 乱 , 将 系数 LH2 和 HL2 利用 序列 1 进行 扩散 加 
密 ， 得 到 加 密 后 的 系数 ;最 后 将 加 密 后 的 系数 矩阵 与 其 他 未 
改变 系数 一 起 进行 小 波 逆 变换 得 到 加 密 后 的 医学 图 像 。 解 密 
过 程 是 加 密 过 程 的 逆 过 程 ， 只 需 将 原始 医学 图 像 改 为 加 密 医 
学 图 像 ， 系 数 LL; 先进 行 混沌 反 置 乱 再 进行 混沌 扩散 ， 其 他 
的 步骤 保持 不 变 即 可 实现 解密 。 算 法 在 公开 的 网 络 环境 下 构 
建 了 一 个 安全 的 图 像 传 输 模型 ， 图 像 的 发 送 者 和 接收 者 协商 
密 钥 ， 并 通过 公 钥 密码 体制 进行 安全 传输 ， 由 于 混沌 加 密 体 
制 具 有 大 的 密 钥 空间 和 强 的 初 值 敏感 性 ， 破 译 者 在 不 知道 密 
钥 情 况 下 仅 从 密 文 中 无 法 破译 出 正确 的 图 像 信息 。 本 文 的 创 
新 之 处 将 2D sine logistic 混沌 映射 应 用 到 了 医学 图 像 的 小 ; 
域 加 密 中 , 并 设计 了 新 的 混沌 扩散 和 置 乱 方法 提高 加 密 效率 ， 
同时 选择 对 LH2 和 HL2 系数 进行 置 乱 加 密 ， 有 效 解决 了 仅 对 
低频 系数 加 密 


邓 小 鸿 ， 


(a) ct (b) 文献 [13 加 罕 结 


图 2 ”加密 效果 对 比 结果 


Fig.2 Comparison results of encrypted efficiency 
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L 一 | 
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吕 混沌 扩散 


> 混沌 置 乱 | > 小 波 逆 变 换 


$n | 
序列 2 
混沌 初 Pe re Pe 加 密 医 
成 


值 /参数 学 图 像 
y 序 列 1 
中 高 频 系数 a 
。 LH;, HL 吕 混沌 扩散 
图 3 本 文 算法 模型 


Fig.3 Model of the proposed algorithm 
2.2 基于 2D sine logistic 的 扩散 和 置 乱 算法 
1) 扩散 机 制 
设 seq, 和 seq, 分 别 为 2D sine logistic 混沌 映射 产生 的 混沌 
序列 ， 其 中 sea, 用 于 对 系数 进行 扩散 ， 而 seq; 用 于 对 扩散 后 的 
位 置 置 乱 。 设 ce 为 原始 医学 图 像 整数 小 波 分 解 后 的 
么 采用 式 (4) 生 成 混沌 密 钥 , 采用 式 (5) 前 向 反馈 
机 制 对 系数 进行 扩散 。 


key(i) = mod(round(seq,(i)*10”), 256) i e[0, len] (4) 


e€_coe(li)= mod((coe(i) @ key(i)) + key(i), 2506) if (i = 

©_coe(i) = mod((coe(i) @ key(i)) + key(i—1),256) Be _coe(li—1) else G) 
其 中 : hey 为 密 钥 ，。-cee 为 加 密 后 系数 集合 ，len 为 系数 集合 
中 元 素 的 个 数 。mod 为 取 余 函数 , round 为 四 舍 五 入 函数 ，@ 
为 异 或 操作 。 
2) 置 乱 机 制 


对 于 扩散 后 的 系数 采用 式 (6) 进 行 置 乱 ， 改 变 系 数位 置 。 
ew _ segq,,location] = sort(segq,) 


@_coe'i)=e_ coe(location(i)) iel[0,M *N] (6) 


其 中 : sorr_segq, 是 对 序列 seq, 进行 升序 排序 后 的 结果 , location 
是 排序 后 元 素 在 原先 序列 中 的 位 置 ，sort 为 升序 排序 函数 ， 
e-coe' 为 对 -cee 进行 位 置 置 乱 后 的 结果 。 
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基于 2D sine logistic 的 扩散 和 置 乱 算法 如 算法 1 所 示 。 
算法 1 根据 生成 的 混沌 序列 对 输入 的 系数 进行 扩散 
和 置 乱 加 密 
输入 : 混沌 序列 集合 seq, 和 segq, ,系数 集合 coe 
输出 : 加 密 系数 集合 e_coe' 
. Ce_coe' =coe; /初始 化 为 原 系数 集合 
.for i=1:len do Wien 为 集合 长 度 
seqi (i)=mod(round( segq, (i)*105),256); 


1 
2 
3 
4. end for 
S$. for i=1:len do 
6. use the equation 5 to encrypt; /扩散 加 密 
7. end for 
8. [sort_seq2,loc]=sort(seq2);// 对 seqz 升序 排序 
9. for i=1:len do 

10.， _ e_coe' (i)=e-_coe (loc(i));// 位 置 置 乱 

11. end for 

12. output e_coe' 

在 算法 1 基础 上 ， 本 文 提出 的 医学 图 像 频 域 加 密 算 法 如 
算法 2 所 示 。 
算法 2 ”根据 给 定 的 混沌 参 

密 医 学 图 像 

输入 : 混沌 初 值 x0,y0, 控 制 参 数 & 和 B ,原始 医学 图 像 I 

输出 加密 医 学 图 像 到 
1. EI1=1; // 初 始 化 为 原始 医学 图 像 
2. [dim1, dim2]=size(D); /获取 医学 图 像 维 数 
3. Level=3;// 设 置 小 波 分 解 的 级 数 
4. decompose53(7, dim1, level);// 小 波 正 变换 
5. N=dim1/(21evel); 
6 
7 
8 


t 


数 和 原始 医学 图 像 得 到 加 


. cof_LL=s33(1:N,1:N);//s33 为 三 级 分 解 的 低频 系数 
. cof_HL=s22(1:N*2, N*2+1:2*N*2);//s22 为 二 级 分 级 的 低频 系数 
. cof_LH=d22(1:N*2, 1:2*N);//d22 为 二 级 分 解 的 高 频 系数 
[seqi,seq2]=D2_SLMM( & , 0 ,x0,y0, 4*N*N);//D2_SLMM 为 用 式 (3) 
生成 混沌 序列 的 函数 
10. for ae cof_LL do 
11. ”调用 算法 1 得 到 e_cof_LL; 
12. end for 
13. for be cof_HL and c ecof_LH do 
14. ”调用 算法 1 中 置 乱 方 法 得 到 e_cof_HL 和 ee_cof_LH:; 
15. end for 
16. s33(1:N,1:N)= e_cof_LL.: 
17. s22(1:N*2, N*2+1:2*N*2)= e_cof_HL; 
18. d22(1:N*2, 1:2*N)= e_cof_LH:;: 
19. EI=recompose53(T’, dim1, level);// 小 波 道 变换 
20. output E/; 


3 ”实验 结果 与 分 析 


实验 中 选取 的 6 幅 图 像 医学 图 像 均 来 自 于 中 南大 学 湘 雅 
医学 院 ， 其 中 4 幅 尺 寸 为 512*$12(x-ray、US 、CT、MRID，2 
幅 尺 寸 为 1024*1024(CT_foot 和 MRI_cervices)。 实 验 主 机 配 
置 环境 : MATLAB 7.0，CPU 为 Intel® CoreTM i5-6500 3.20 
GHz， 内 存 8 GB，64 位 Windows 7 旗舰 版 操作 系统 。 实 验 
中 设 定 2D sine Logistic 混沌 的 初 值 x0=0.9380，y0=0.7006， 
控制 参数 “=1,2=3。 

3.1 加 密 结果 及 分 析 
3.1.1 加 密 结果 
采用 本 文 算法 对 载体 医学 图 像 进行 加 密 后 的 结果 如 


党 


路 
上 
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所 示 。 图 


邓 小 鸿 ， 等 : 基于 2DSineLogistic 混沌 映射 的 医学 图 像 频 域 加 密 算 法 


4 中 (ac)(e)(g)G)(o 分 别 为 原始 医学 图 


人 b)(d)(GDG)G)G) 为 其 对 应 的 加 密 医 学 图 像 。 从 图 4 中 所 示 结 果 


来 看 ， 加 密 图 像 具 


随机 噪声 分 布 特 


像 ， 


性 ， 具 有 较 好 的 力 


果 , 当 加 密 图 像 未 受到 窜改 和 解密 密 钥 与 加 密 时 保持 一 


三 


能 无 损 解 密 出 原始 


医学 图 像 。 另 外 ， 


对 比 图 


结果 ， 可 以 明显 发 


问题 。 本 文 算法 通 


1 密 效 
致 时 ， 


2(b) 和 图 4(d) 中 
山本 文 算法 有 效 解决 了 加 密 图 像 中 的 轮廓 


过 分 析 小 波 分 解 特点 ， 选 取 二 级 小 波 分 级 


后 的 中 高 频 子 带 (代表 了 图 像 的 纹理 和 边缘 信息 ) 进行 加 密 ， 


这 样 加 密 后 的 系数 在 进行 小 波 逆 变 换 时 进 


构 ， 充 分 掩盖 图 像 


的 边缘 信息 。 


行 系数 迭 代 和 重 
实验 中 同样 对 三 级 分 解 的 中 


高 频 子 带 进 行 加 密 ， 发 现 效 果 不 如 二 级 系数 明显 ， 选 用 一 级 


分 解 的 中 高 频 系 数 进行 加 密 ， 虽 然 能 较 好 地 解决 轮廓 问题 ， 
但 加 密 时 间 显 著 增 加 。 
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医学 图 像 加 密 结 果 

Fig.4 Encryption results of medical images 
3.1.2 性 能 分 析 
1) 密 钥 安 全 分 析 
密 钥 安全 性 主要 通过 密 钥 空间 和 密 钥 敏 感性 来 衡量 ， 密 
钥 空间 大 小 决定 了 攻击 者 采用 穷 举 法 进行 暴力 破解 的 难度 ， 
而 密 钥 敏感 性 确保 破译 者 不 能 使 用 与 实际 密 钥 相近 的 密 钥 恢 
复出 原始 信息 。 首 先 本 文 设计 的 密 钥 长 度 为 128 位 (混沌 初 值 
和 控制 参数 分 别 占 4 个 字 节 ， 密 钥 位 数 为 4*32=128 位 )， 则 
密 钥 空间 为 23”， 通 过 穷 举 搜索 法 破解 密 钥 从 计算 上 是 不 可 
能 的 。 其 次 通过 实验 测试 ， 将 2D sine logistic 混沌 的 初 值 x0 
修改 为 0.9381， 其 他 的 参数 不 变 ， 采 用 该 密 钥 去 解密 图 像 ， 
得 到 的 测试 结果 如 图 5 所 示 。 图 5 中 ，(a) 为 原始 医学 图 像 ， 
(b) 为 采用 默认 密 钥 进行 加 密 后 的 图 像 , (c) 为 采用 修改 密 钥 进 
行 加 密 后 的 图 像 ,， (qd) 为 (b) 和 (c) 中 图 像 求 差 的 结果 。 从 图 5(d) 
中 可 以 看 出 ， 密 钥 的 微小 变化 会 造成 加 密 结果 的 显著 不 同 ， 
另外 采用 修改 后 的 密 钥 进行 解密 也 无 法 得 到 原始 的 医学 图 像 。 
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图 5 密 钥 敏感 性 结果 

Fig.5 The results of key sensitivity 
2) 抗 统计 攻击 分 析 
加 密 图 像 信息 炉 和 相 邻 像素 相关 性 是 定量 衡量 算法 抗 统 


计 攻 击 的 有 效 指标 。 图 像 信 息 炉 表示 图 像 灰 度 分 布 的 统计 特 
征 ， 对 于 一 幅 类 似 随 机 噪声 的 加 密 图 像 ， 其 灰 度 分 布 应 该 趋 
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向 于 均匀 化 ， 其 期 望 值 应 该 为 8。 信 息 焙 和 像素 相关 性 计算 


公式 如 下 : 


Entropy = -> pll) log (p(D) 


EX 一 人 他 一 0 


oo 
x 


(7) 


(8) 


式 (7) 中 ，N 为 像素 值 的 灰 度 级 数 ，p( 为 像素 值 ' 出现 的 


概率 。 式 (8) 中 ,EE 为 求 期 望 , X 和 了 为 两 组 数据 序列 ， 


Hx 和 


必 分 别 为 数据 序列 的 均值 ，ox 和 为 标准 差 。 如 果 T 
列 X 和 节 具 有 高 的 相关 性 ，corr 
经 过 计算 ， 加 


个 序 
值 接近 于 1， 相 反 接近 于 0。 
密 后 的 6 幅 图 像 的 信息 烂 在 7.8 左右 ， 接 近 于 
理想 值 8。 表 2 给 出 了 载体 图 像 CT_foot 加 密 前 后 像素 的 相 


关 性 。 从 表 2 中 结果 可 以 看 出 ， 加 密 图 像 具 有 非常 低 的 像素 


相关 性 。 综 上 ， 本 文 加 密 算法 能 较 好 地 抵抗 统计 攻击 。 
表 2 像素 相关 性 测试 结果 
Table 2 Testing results of pixel correlation 
载体 图 像 水 平方 向 垂直 方向 对 角 线 方向 
CT_foot 原始 图 像 0.9988 0.9661 0.95942 
CT _foot 加 密 图 像 0.0017 -0.0034 0.0308 


3) 抗 差分 攻击 分 析 
差分 攻击 是 通过 比较 分 析 


特定 区 别 的 明文 在 通过 加 密 


后 的 变化 传播 情况 来 攻击 密码 算法 。 加 密 图 像 对 明文 的 敏感 


性 强度 是 衡量 算法 抗 差分 攻击 的 常 
中 的 像素 改变 率 NPCR 计算 。 


和 
NPCR = D2 DG. x100% 
MN 


DEND=CCD=C(0,)7:1,0) 


| 方法， 通常 利用 


式 (9) 


(9) 


其 中 : ww*w 为 图 像 的 尺寸 ，“ 和 < 分 别 是 两 幅 医 学 图 像 石和 
己 对 应 的 密 文 ， 其 中 五 和 五 仅 有 一 个 像素 值 不 同 。 

0.9975 
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0.9965 
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0.9955 


0.995 
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6 ”像素 改变 率 测 试 结果 
Fig.6 Testing results of pixel changing rate 


6 给 出 了 NPCR 值 测试 结果 。 测 试 中 选 


] 60 组 医学 


图 像 (每 副 载体 图 像 修 改 10 个 不 同位 置 的 像素 点 ， 与 原 
成 一 组 )， 对 每 组 进行 编号 ， 并 随机 选择 每 组 进行 加 密 ， 


想 值 0.99603 上 下 波动 ， 说 明 算 法 具有 较 强 的 抗 差分 攻 
Ne 
4 加 解密 效率 分 析 


表 3 给 


看 出 ， 当 医学 图 像 尺 寸 为 512*512 时 ，4 幅 载 体 图 像 的 
1 


计算 每 组 的 VPCR 值 。 从 图 6 中 可 以 看 出 ， 像 素 改变 率 在 理 


出 了 6 副 载 体 图 像 加 密 和 解密 所 需 的 时 间 ， 表 中 
数据 均 是 对 图 像 加 解密 10 次 得 到 的 平均 结果 。 从 表 3 中 可 以 
学 加 解 


图 组 
分 别 


ee 
证 用 


密 时 间 在 1.1 s 左右 ， 当 尺寸 为 1024*1024 时 ，2 幅 载体 
的 加 解密 时 间 在 5s 左右 ， 算 法 加 解密 效率 较 好 。 
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表 3 本文 算法 加 解密 效率 测试 结果 以 上 两 个 文献 进行 对 比 ， 对 比 的 性 能 指标 包括 加 密 图 像 信息 
Table 3 Encrypt efficiency testing results of the proposed algorithm 粮 、 加 密 时 间 和 加 密 图 像 是 否 有 轮 廊 。 表 4 给 出 了 本 文 算法 
载体 图 像 加 密 时 间 /s 解密 时 间 /s 和 文献 [13,15] 的 对 比 结果 。 从 表 4 中 结果 可 以 看 出 ， 本 文 算 
x-ray(512*512) 1.129 1.131 法 相 比 文献 [15] 方 法 ， 虽 然 加 密 图 像 信息 焙 较 低 ， 但 在 加 密 
us(512*512) 1.137 1.132 时 间 上 具有 明显 优势 ， 对 于 一 副 M*N 的 图 像 ， 文献 [15] 需 要 
ct(512*512) 1.13 1.129 加 密 的 像素 个 数 为 M*N, 而 本 文 算 法 需要 加 密 的 系数 个 数 仅 
mri(512*512) 1.14 1.396 为 M*N(1/256+1/23) 守 M*N/8。 与 文献 [13] 相 比 ， 本 文 算法 的 加 
CT_foot(1024*1024) 4.966 4.847 密 图 像 信息 粒 更 接近 于 理想 值 8， 因 为 选取 了 更 多 的 系数 进 
MRI_cervices(1024*1024) 4.878 4.952 行 混沌 加 密 ， 但 加 密 时 间 上 两 个 算法 差不多 ， 原 因 在 于 文献 
3.2 算法 对 比 结果 及 分 析 [13] 方 法 生成 超 混沌 序列 的 复杂 度 高 于 2D sine logistic。 除 了 
本 文 算法 思路 主要 来 源 于 文献 [13,15]， 实 验 结果 主要 与 言 息 灶 外 ， 本 文 算法 最 大 的 优势 是 解决 了 加 密 图 像 中 的 轮 廊 问 题 。 

表 4 本 文 算法 与 其 他 对 比方 法 的 测试 结果 


Table 4 Comparison results between the proposed algorithm and other schemes 


载体 图 像 文献 [15] 文献 [13] 本 文 算法 
加 密 图 像 信 ”加密 图 像 信 ，,，，， 密 图 像 信 ，,， ，， 是 否 有 轮 

加 密 时 间 s 是 否 有 轮 廊 加 密 时 间 。 是 否 有 轮廓 加 密 时 间 s 

娠 炉 乱 业 姑 入 廊 

X-ray 7.998 41.256 无 7.792 1.096 无 7.803 1.129 无 

us 7.997 40.977 无 7.794 0.997 无 7.812 1.137 无 

ct 7.998 41.023 无 7.801 1.002 有 7.810 1.131 无 

mri 7.996 41.144 无 7.800 1.078 有 7.825 1.145 无 

CT_foot 7.999 198.56 无 7.796 5.002 有 7.799 4.966 无 

MRI cervices 7.998 199.87 无 7.795 4.993 有 7.803 4.878 无 
§ for medical applications [J]. Communications in Nonlinear Science and 

4 ”结束 语 


Numerical Simulations, 2015, 24(1): 98-116. 


网 络 环境 下 医学 图 像 的 安全 传输 需求 日 益 增加 ， 基 于 泥 
沌 的 医学 图 像 加 密 方法 受到 研究 人 员 的 广泛 关注 。 本 文 在 分 
析 现 有 方法 的 基础 上 ， 充 分 考虑 医学 图 像 自 身 特点 ， 提 出 了 
基于 2D sine Logistic 混沌 映射 的 医学 图 像 频 域 加 密 算法 ， 设 
计 了 新 的 混沌 扩散 和 置 乱 算法 ， 对 三 级 小 波 分 解 的 低频 系数 
行 加 密 ， 提 高 了 算法 的 安全 性 和 效率 ， 另 外 为 了 解决 加 密 
图 像 中 出 现 的 轮廓 问题 ， 选 取代 表 边 缘 和 纹理 信息 的 二 级 小 
波 分 级 的 中 高 频 系数 进行 混沌 扩散 。 实 验 结果 表明 本 文 算 法 
在 加 密 安 全 性 和 效率 方面 具有 优势 ， 算 法 适用 于 大 数据 量 的 
医学 图 像 实时 加 密 。 
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